Combined Spectroscopic, Thermodynamic, and Theoretical Approach for Detecting and Quantifying Hydrogen Bonding and Dispersion Interaction in Ionic Liquids.

ConspectusIonic liquids (ILs) are attracting increasing interest in science and engineering due to their unique properties that can be tailored for specific applications. Clearly, a better understanding of their behavior on the microscopic scale will help to elucidate macroscopic fluid phenomena and thereby promote potential applications. The advantageous properties of these innovative fluids arise from the delicate balance of Coulomb interactions, hydrogen bonding, and dispersion forces. The development of these properties requires a fundamental understanding of the strength, location, and direction of the different types of interactions and their contribution to the overall phase behavior. Contrary to expectations, hydrogen bonding and dispersion interactions have a significant influence on the structure, dynamics, and phase behavior of ILs.The synergy between experimental and theoretical methods has now advanced to a stage where hydrogen bonds and dispersion effects as well as the competition between the two can be studied in detail. In this account, we demonstrate that a suitable combination of spectroscopic, thermodynamic, and theoretical methods enables the detection, dissection, and quantification of noncovalent interactions, even in complex systems such as ionic liquids. This approach encompasses far-infrared vibrational spectroscopy (FIR), various thermodynamic methods for determining enthalpies of vaporization, and quantum chemical techniques that allow us to switch dispersion contributions on or off when calculating the energies and spectroscopic properties of clusters.We briefly discuss these experimental and theoretical methods, before providing various examples illustrating how the mélange of Coulomb interaction, hydrogen bonds, and dispersion forces can be analyzed, and their individual contributions quantified. First, we demonstrated that both hydrogen bonding and dispersion interactions are manifested in the FIR spectra and can be quantified by observed shifts of characteristic spectral signatures. Through the selection of suitable protic ionic liquids (PILs) featuring anions with varying interaction strengths and alkyl chain lengths, we were able to demonstrate that dispersion interactions can compete with hydrogen bonding. The resultant transition enthalpy serves as a measure of the dispersion interaction. Contrary to expectations, PILs possess lower enthalpies of vaporization compared with aprotic ILs (AILs). The reason for this is simple: In protic ILs, ion pairs carry both the hydrogen bond and attractive dispersion between the cation and anion into the gas phase. By utilizing a well-curated set of protic ILs and molecular analogues, we successfully disentangled Coulomb interaction, hydrogen bonding, and dispersion interaction through purely thermodynamic methods.

Role of Hydrogen Bond Defects for Cluster Formation and Distribution in Ionic Liquids by Means of Neutron Diffraction and Molecular Dynamics Simulations.

Defects fundamentally govern the properties of all real materials. Correlating molecular defects to macroscopic quantities remains a challenge, particularly in the liquid phase. Herein, we report the influence of hydrogen bonds (HB) acting as defects in mixtures of non-hydroxyl-functionalized ionic liquids (ILs) with an increasing concentration of hydroxyl-functionalized ILs. We observed two types of HB defects: The conventional HBs between cation and anion (c-a), and the elusive HBs between cations (c-c) despite the repulsive Coulomb forces. We use neutron diffraction with isotopic substitution in combination with molecular dynamics simulations for measuring the geometry, strength, and distribution of mobile OH defects in the IL mixtures. In principle, this procedure allows relating the number and stability of defects to macroscopic properties such as diffusion, viscosity, and conductivity, which are of utmost importance for the performance of electrolytes in batteries and other electrical devices.

The Influence of Deuterium Isotope Effects on Structural Rearrangements, Ensemble Equilibria, and Hydrogen Bonding in Protic Ionic Liquids.

We report strong isotope effects for the protic ionic liquid triethylammonium methanesulfonate [TEA][OMs] by means of deuterium solid-state NMR spectroscopy covering broad temperature ranges from 65 K to 313 K. Both isotopically labelled PILs differ in non-deuterated and fully deuterated ethyl groups of the triethyl ammonium cations. The N-D bond of both cations is used as sensitive probe for hydrogen bonding and structural ordering. The 2 H NMR line shape analysis provides the deuteron quadrupole coupling constants and the characteristics of a broad heterogeneous phase with simultaneously present static and mobile states indicating plastic crystal behavior. The temperatures where both states are equally populated differ by about 80 K for the two PILs, showing that deuteration of the ethyl groups in the trialkylammonium cations tremendously shifts the equilibrium towards the static state. In addition, it leads to a significant less cooperative transition, associated with a significantly reduced standard molar transition entropy.

Structure, hydrogen bond dynamics and phase transition in a model ionic liquid electrolyte.

We show that solid-state NMR spectroscopy is a suitable method for characterizing the structure, hydrogen bond dynamics and phase transition behavior in protic ionic liquids (PILs). Deuteron line shape and spin relaxation time analysis provide a description of the structural and dynamical heterogeneity in the solid state of the model PIL triethyl ammonium bis(trifluoromethanesulfonyl)amide [TEA][NTf2]. Therein, we observed two deuteron quadrupole coupling constant for the ND bond of the TEA cation, indicating differently strong hydrogen bonds to the nitrogen and oxygen atoms of the NTf2 anion, as we could confirm by DFT calculations. The transition processes in the dynamically heterogeneous phase are characterized by two standard molar enthalpies and thus different stages of melting. We provide geometry, rates and energetics of the cation in the solid and liquid states of the PIL. Comparison with PILs having stronger interacting anions shows higher enthalpy change between the solid and liquid states, lower activation barriers of tumbling motion and higher amplitude of librational motion for the TEA cation in the presence of the weakly interacting anion NTf2. We provide reasonable relations between microscopic and macroscopic properties, as is relevant for any kind of application.

Three in One: The Versatility of Hydrogen Bonding Interaction in Halide Salts with Hydroxy-Functionalized Pyridinium Cations.

The paradigm of supramolecular chemistry relies on the delicate balance of noncovalent forces. Here we present a systematic approach for controlling the structural versatility of halide salts by the nature of hydrogen bonding interactions. We synthesized halide salts with hydroxy-functionalized pyridinium cations [HOCn Py]+ (n=2, 3, 4) and chloride, bromide and iodide anions, which are typically used as precursor material for synthesizing ionic liquids by anion metathesis reaction. The X-ray structures of these omnium halides show two types of hydrogen bonding: 'intra-ionic' H-bonds, wherein the anion interacts with the hydroxy group and the positively charged ring at the same cation, and 'inter-ionic' H-bonds, wherein the anion also interacts with the hydroxy group and the ring system but of different cations. We show that hydrogen bonding is controllable by the length of the hydroxyalkyl chain and the interaction strength of the anion. Some molten halide salts exhibit a third type of hydrogen bonding. IR spectra reveal elusive H-bonds between the OH groups of cations, showing interaction between ions of like charge. They are formed despite the repulsive interaction between the like-charged ions and compete with the favored cation-anion H-bonds. All types of H-bonding are analyzed by quantum chemical methods and the natural bond orbital approach, emphasizing the importance of charge transfer in these interactions. For simple omnium salts, we evidenced three distinct types of hydrogen bonds: Three in one!

Balance Between Contact and Solvent-Separated Ion Pairs in Mixtures of the Protic Ionic Liquid [Et3NH][MeSO3] with Water Controlled by Water Content and Temperature.

The formation of aggregates of ionic species is a crucial process in liquids and solutions. Ion speciation is particularly interesting for the case of ionic liquids (ILs) since these Coulombic fluids consist solely of ions. Most of their unique properties, such as enthalpies of vaporization and conductivities, are strongly related to ion pair formation. Here, we show that the balance of hydrogen-bonded contact ion pairs (CIP) and solvent-separated (SIP) ion pairs in protic ionic liquids (PILs) and in their mixtures with water can be well understood by a combination of far-infrared (FIR) and mid-infrared (MIR) spectroscopy, density functional theory (DFT) calculations of PIL/water aggregates, and molecular dynamics (MD) simulations of PIL/water mixtures. This combined approach is applied to mixtures of triethylammonium methanesulfonate [Et3NH][MeSO3] with water. It is shown that ion speciation in this mixture depends on three parameters: the relative hydrogen bond acceptor strength of the counter ion and the molecular solvent, the solvent concentration, and the temperature. For selected PIL/water mixtures, the equilibrium constants for CIPs and SIPs were determined as a function of the solvent content and temperature. Finally, for the studied PIL/water mixtures, the transition from CIPs to SIPs could be understood on enthalpic and entropic grounds. A detailed picture of this interconversion process could be described at the molecular level by means of MD simulations. In addition, the concentration dependence of ion pair formation can be well understood with help of a simplified "cartoon-like" statistical model describing hydrogen bond redistribution.

Freezing the Motion in Hydroxy-Functionalized Ionic Liquids-Temperature Dependent NMR Deuteron Quadrupole Coupling Constants for Two Types of Hydrogen Bonds Far below the Glass Transition.

We measured the deuteron quadrupole coupling constants (DQCCs) for hydroxy-functionalized ionic liquids (ILs) with varying alkyl chain length over the temperature range between 60 and 200 K by means of solid-state NMR spectroscopy. For all temperatures, the 2H spectra show two DQCCs representing different types of hydrogen bonds. Higher values, ranging from 220 to 250 kHz, indicate weaker hydrogen bonds between cation and anion (c-a), and lower values varying from 165 to 210 kHz result from stronger hydrogen bonds between the OD groups of cations (c-c), in agreement with recent observations in infrared, neutron diffraction, and NMR studies. We observed different temperature dependencies for (c-a) and (c-c) hydrogen bonding. From the static pattern of the 2H spectra at the lowest temperatures, we derived the true DQCCs being up to 20 kHz larger than recently reported values measured at the glass transition temperature. We were able to freeze the librational motions of the hydrogen bonds in the ILs. The temperature dependence of the (c-a) and (c-c) cluster populations in the glassy state is opposite to that observed in the liquid state, partly anticipating the behavior of ILs tending to crystallize.

Counting cations involved in cationic clusters of hydroxy-functionalized ionic liquids by means of infrared and solid-state NMR spectroscopy.

In hydroxy-functionalized ionic liquids, two types of hydrogen bonding coexist: the conventional H-bonds between cation and anion (c-a) and those between cation and cation (c-c), although the interaction between like-charged ions is supposed to be much weaker due to the repulsive Coulomb forces. Counting the cations involved in either (c-a) or (c-c) clusters is a challenge. For that purpose, we recently performed neutron diffraction (ND) measurements and molecular dynamics (MD) simulations at and above room temperature accompanied by NMR solid-state experiments in the glassy state of the ILs. In principle, these methods are suitable for determining the populations of (c-a) and (c-c) cluster species. For different reasons we could only address single temperatures and/or small temperature intervals above 300 K. The by far largest temperature range with reasonable efforts is accessible by simple infrared (IR) spectroscopy. However, counting (c-a) or (c-c) hydrogen bonds is a difficult task due to the different transition dipole moments resulting in varying intensities and broad vibrational bands. Here we present a method for deriving the number of cations involved in (c-a) ion pairs from IR spectra in the OH stretch region. This procedure provides access to the equilibria of (c-a) and (c-c) hydrogen bonds as a function of temperature allowing derivation of the transition enthalpy.

Controlling "like-likes-like" charge attraction in hydroxy-functionalized ionic liquids by polarizability of the cations, interaction strength of the anions and varying alkyl chain length.

We provide comprehensive understanding of "like-likes-like" charge attraction in hydroxy-functionalized ionic liquids (ILs) by means of infrared spectroscopy (IR), quantum chemistry and differential scanning calorimetry (DSC). We show that hydrogen bonding between cation and cation (c-c) is possible despite the repulsive forces between ions of like charge. Already at room temperature, the (c-c) hydrogen bonds can compete with the regular Coulomb-enhanced hydrogen bonds between cation and anion (c-a). For a large set of well-selected ILs, we show that "like-charge attraction" between the OH-functionalized cations is controllable by the polarizability of the cation, the interaction strength of the anion and the length of the hydroxyalkyl chain. In particular, we clarify whether tethering the OH group away from the positive charge center of the cationic ring with longer hydroxyalkyl chains compensates for unfavourable cation/anion combinations with respect to (c-c) cluster formation. For that purpose, we synthesized and characterized twelve ionic liquids including the differently polarizable cations, 1-(n-hydroxyalkyl)-1-methylpiperidinium [HOCnMPip]+ and 1-(n-hydroxyalkyl)-pyridinium [HOCnPy]+, as well as the weakly and strongly interacting anions, bis(trifluoromethanesulfonyl)imide [NTf2]- and methanesulfonate [OMs]-, respectively. On top, we varied the hydroxyalkyl chain length (HOCn) (n = 2-5). We systematically show how these three molecular ion parameters affect like-charge attraction. The use of polarizable cations, weakly interacting anions, and long alkyl chain tethers results in (c-c) clustering already at room temperature. Kinetic trapping is not a prerequisite for the existence of (c-c) cluster species in ILs. Moreover, we demonstrate that micro structuring affects macroscopic behavior of this type of ILs. We observed that substantial (c-c) interaction prevents ILs from crystallizing. Instead, these ILs supercool and finally form a glass.

Hydrogen Bonding Between Ions of Like Charge in Ionic Liquids Characterized by NMR Deuteron Quadrupole Coupling Constants-Comparison with Salt Bridges and Molecular Systems.

We present deuteron quadrupole coupling constants (DQCC) for hydroxyl-functionalized ionic liquids (ILs) in the crystalline or glassy states characterizing two types of hydrogen bonding: The regular Coulomb-enhanced hydrogen bonds between cation and anion (c-a), and the unusual hydrogen bonds between cation and cation (c-c), which are present despite repulsive Coulomb forces. We measure these sensitive probes of hydrogen bonding by means of solid-state NMR spectroscopy. The DQCCs of (c-a) ion pairs and (c-c) H-bonds are compared to those of salt bridges in supramolecular complexes and those present in molecular liquids. At low temperatures, the (c-c) species successfully compete with the (c-a) ion pairs and dominate the cluster populations. Equilibrium constants obtained from molecular-dynamics (MD) simulations show van't Hoff behavior with small transition enthalpies between the differently H-bonded species. We show that cationic-cluster formation prevents these ILs from crystallizing. With cooling, the (c-c) hydrogen bonds persist, resulting in supercooling and glass formation.

The Double-Faced Nature of Hydrogen Bonding in Hydroxy-Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations.

We characterize the double-faced nature of hydrogen bonding in hydroxy-functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H⋅⋅⋅O and O⋅⋅⋅O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation-cation interaction is stronger than the cation-anion interaction. We compare the hydrogen-bond geometries of both "doubly charged hydrogen bonds" with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like-charge attraction is almost competitive with conventional ion-pair formation.

Dissecting the Vaporization Enthalpies of Ionic Liquids by Exclusively Experimental Methods: Coulomb Interaction, Hydrogen Bonding, and Dispersion Forces.

The quantification of hydrogen bonding and dispersion energies from vaporization enthalpies is a great challenge. Dissecting interaction energies is particularly difficult for ionic liquids (ILs), for which the composition of the different types of interactions is known neither for the liquid nor for the gas phase. In this study, we demonstrate the existence of ion pairs in the gas phase and dissect the interaction energies exclusively from measured vaporization enthalpies of different alkylated protic ILs (PILs) and aprotic ILs (AILs) and the molecular analogues of their cations. We demonstrate that the evaporated ion pairs are characterized by H-bond-enhanced Coulomb interaction. The overall interaction energy for the ILs in the bulk phase is composed of Coulomb interaction (76 kJ mol-1 ), hydrogen bonding (38 kJ mol-1 ), and minor dispersion interaction (10 kJ mol-1 ). Thus, hydrogen bonding prominently contributes to the overall interaction energy of PILs, which is reflected in the properties of this class of liquids.

Like-likes-Like: Cooperative Hydrogen Bonding Overcomes Coulomb Repulsion in Cationic Clusters with Net Charges up to Q=+6e.

Quantum chemical calculations have been employed to study kinetically stable cationic clusters, wherein the monovalent cations are trapped by hydrogen bonding despite strongly repulsive electrostatic forces. We calculated linear and cyclic clusters of the hydroxy-functionalized cation N-(3-hydroxypropyl) pyridinium, commonly used as cation in ionic liquids. The largest kinetically stable cluster was a cyclic hexamer that very much resembles the structural motifs of molecular clusters, as known for water and alcohols. Surprisingly, strong cooperative hydrogen bonds overcome electrostatic repulsion and result in cationic clusters with a high net charge up to Q=+6e. The structural, spectroscopic, and electronic signatures of the cationic and related molecular clusters of 3-phenyl-1-propanol could be correlated to NBO parameters, supporting the existence of "anti-electrostatic" hydrogen bonds (AEHB), as recently suggested by Weinhold. We also showed that dispersion forces enhance the cationic cluster formation and compensate the electrostatic repulsion of one additional positive charge.

Dynamical heterogeneities in ionic liquids as revealed from deuteron NMR.

The heterogeneity in dynamics has important consequences for understanding the viscosity, diffusion, ionic mobility, and the rates of chemical reactions in technology relevant systems such as polymers, metallic glasses, aqueous solutions, and inorganic materials. Herein, we study the spatial and dynamic heterogeneities in ionic liquids by means of solid state NMR spectroscopy. In the 2H spectra of the protic ionic liquid [TEA][OTf] we observe anisotropic and isotropic signals at the same time. The spectra measured below the melting temperature at 306 K could be simulated by a superposition of the solid and liquid line shapes, which provided the transition enthalpies between the rigid and mobile fractions. Consequently, we measured the spin-lattice relaxation times T1 for the anisotropic and the isotropic signals for the temperature range between 203 and 436 K. Both dispersion curves could be fitted to models including rotational correlation times, activation barriers and rate constants. This approach allowed determining the rotational correlation times for the N-D molecular vector of the [TEA]+ cation in differently mobile environments. The mobility is only slightly different, as indicated by small differences in activation energies for these processes. The NMR correlation times for the highly mobile phase are linearly related to measured viscosities, which supports the applicability of the Stokes-Einstein-Debye relation.

Characterization of Doubly Ionic Hydrogen Bonds in Protic Ionic Liquids by NMR Deuteron Quadrupole Coupling Constants: Differences to H-bonds in Amides, Peptides, and Proteins.

We present the first deuteron quadrupole coupling constants (DQCCs) for selected protic ionic liquids (PILs) measured by solid-state NMR spectroscopy. The experimental data are supported by dispersion-corrected density functional theory (DFT-D3) calculations and molecular dynamics (MD) simulations. The DQCCs of the N-D bond in the triethylammonium cations are the lowest reported for deuterons in PILs, indicating strong hydrogen bonds between ions. The NMR coupling parameters are compared to those in amides, peptides, and proteins. The DQCCs show characteristic behavior with increasing interaction strength of the counterion and variation of the H-bond motifs. We report the similar presence of the quadrupolar splitting pattern and the narrow liquid line in the NMR spectra over large temperature ranges, indicating the heterogeneous nature of PILs.

Dispersion and Hydrogen Bonding Rule: Why the Vaporization Enthalpies of Aprotic Ionic Liquids Are Significantly Larger than those of Protic Ionic liquids.

It is well known that gas-phase experiments and computational methods point to the dominance of dispersion forces in the molecular association of hydrocarbons. Estimates or even quantification of these weak forces are complicated due to solvent effects in solution. The dissection of interaction energies and quantification of dispersion interactions is particularly challenging for polar systems such as ionic liquids (ILs) which are characterized by a subtle balance between Coulomb interactions, hydrogen bonding, and dispersion forces. Here, we have used vaporization enthalpies, far-infrared spectroscopy, and dispersion-corrected calculations to dissect the interaction energies between cations and anions in aprotic (AILs), and protic (PILs) ionic liquids. It was found that the higher total interaction energy in PILs results from the strong and directional hydrogen bonds between cation and anion, whereas the larger vaporization enthalpies of AILs clearly arise from increasing dispersion forces between ion pairs.

Spectroscopic Evidence for Clusters of Like-Charged Ions in Ionic Liquids Stabilized by Cooperative Hydrogen Bonding.

Direct spectroscopic evidence for hydrogen-bonded clusters of like-charged ions is reported for ionic liquids. The measured infrared O-H vibrational bands of the hydroxyethyl groups in the cations can be assigned to the dispersion-corrected DFT calculated frequencies of linear and cyclic clusters. Compensating the like-charge Coulomb repulsion, these cationic clusters can range up to cyclic tetramers resembling molecular clusters of water and alcohols. These ionic clusters are mainly present at low temperature and show strong cooperative effects in hydrogen bonding. DFT-D3 calculations of the pure multiply charged clusters suggest that the attractive hydrogen bonds can compete with repulsive Coulomb forces.

Controlling the subtle energy balance in protic ionic liquids: dispersion forces compete with hydrogen bonds.

The properties of ionic liquids are determined by the energy-balance between Coulomb-interaction, hydrogen-bonding, and dispersion forces. Out of a set of protic ionic liquids (PILs), including trialkylammonium cations and methylsulfonate and triflate anions we could detect the transfer from hydrogen-bonding to dispersion-dominated interaction between cation and anion in the PIL [(C6 H13 )3 NH][CF3 SO3 ]. The characteristic vibrational features for both ion-pair species can be detected and assigned in the far-infrared spectra. Our approach gives direct access to the relative strength of hydrogen-bonding and dispersion forces in a Coulomb-dominated system. Dispersion-corrected density functional theory (DFT) calculations support the experimental findings. The dispersion forces could be quantified to contribute about 2.3 kJ mol(-1) per additional methylene group in the alkyl chains of the ammonium cation.

Ion pairing in protic ionic liquids probed by far-infrared spectroscopy: effects of solvent polarity and temperature.

The cation-anion and cation-solvent interactions in solutions of the protic ionic liquid (PIL) [Et3NH][I] dissolved in solvents of different polarities are studied by means of far infrared vibrational (FIR) spectroscopy and density functional theory (DFT) calculations. The dissociation of contact ion pairs (CIPs) and the resulting formation of solvent-separated ion pairs (SIPs) can be observed and analyzed as a function of solvent concentration, solvent polarity, and temperature. In apolar environments, the CIPs dominate for all solvent concentrations and temperatures. At high concentrations of polar solvents, SIPs are favored over CIPs. For these PIL/solvent mixtures, CIPs are reformed by increasing the temperature due to the reduced polarity of the solvent. Overall, this approach provides equilibrium constants, free energies, enthalpies, and entropies for ion-pair formation in trialkylammonium-containing PILs. These results have important implications for the understanding of solvation chemistry and the reactivity of ionic liquids.

Equilibrium of contact and solvent-separated ion pairs in mixtures of protic ionic liquids and molecular solvents controlled by polarity.

Polarity controls the equilibrium constants and free energies of contact ion pairs (CIPs) and solvent-separated ion pairs (SIPs) in mixtures of protic ionic liquids and molecular solvents. The subtle balance between the ionic species was studied by far-infrared difference spectra and related DFT-calculated properties for solvents of low and high polarity and for different solvent concentrations.